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Abstract. In this article we explore the origin of black hole thermodynamics using Semi- 
co classical states in loop quantum gravity. We re-examine the case of entropy using a density 
matrix for a coherent state and describe correlations across the horizon due to SU(2) in- 
tertwiners. We further show that Hawking radiation is a consequence of a non-Hermitian 
term in the evolution operator, which is necessary for entropy production or depletion at the 
""TT horizon. This non-unitary evolution is also rooted in formulations of irreversible physics. 
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1 Introduction 
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With new results in lattice gravity, loop quantum gravity (LQG) and various other approaches to 

CN quantum gravity, there is hope that we might find the quantum of space-time [2,4,21,22,39,41]. 

Though unifying gravity with other forces of nature will remain a puzzle. General theory of 

relativity is a theory of gravity and quantum field theory is the theory of particle dynamics 

and both are experimentally verified. Yet when we describe quantum fields in gravitational 

backgrounds, new phenomena arise which seem to signify the existence of new physics. Black 

^^ hole thermodynamics and Hawking radiation are particularly interesting phenomena. The laws 

^vj of black hole mechanics derived in [7] were the first indication that the black hole might be 

T— I a system similar to a thermodynamic ensemble. The area increasing theorem [26] led to the 

^ conjecture that the entropy of the thermodynamic system is proportional to the area of the 

horizon [9] . The area of the horizon divided by four Planck length squared today is known as 

the Bekenstein-Hawking entropy. The use of quantum fields near the horizon then showed that 

black holes emit radiation in thermal spectrum at a temperature given by the surface gravity 

of the black hole [27], and the Planck's constant naturally got included in the quantification of 

entropy, temperature etc. This confirmed that there is a 'quantisation' of the space-time and 

the microscopic structure will explain the origin of thermodynamic quantities. There are several 

explanations of the nature of the microscopic structure of the black hole entropy [4,39], and 

the origin of thermodynamics, some of which will be reported in this volume. In this article 

I shall discuss the use of 'coherent states' in LQG to derive the physics of the system [25,42]. 

The use of coherent states is well known in quantum electrodynamics in particular for lasers, 

and we will use similar coherent states for gravity. Once the coherent states are identified for 

the black hole, the derivation of thermodynamics follows the standard entropy formulation of 

correlated systems, part of which has been traced away [18]. A reduced density matrix is derived 

for a region outside the horizon and the wavefunction inside the horizon is traced over. This 
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reduced density matrix gives an entropy proportional to the area of the horizon. There are 
semiclassical corrections to the Bekenstein-Hawking term, and the nature of the corrections is 
a function of the discretisation or graph embeddings used to describe the system [19]. In this 
article we describe the derivation of entropy using coherent states, and discuss briefly the nature 
of correlations which arise due to the gauge invariant coherent states using SU(2) intertwiners. 
These calculations are a preliminary indication of what a 'physical coherent state' might carry 
as correlations at the horizon [15]. The correlations make the state inside the horizon 'entangled' 
with the state outside the horizon. 

As the coherent state is defined in the canonical formalism where space-time is foliated with 
constant time slices, the coherent state in discussion is defined in one time slice, which we 
take as the initial time slice. We then use a semi-classical Hamiltonian to evolve the coherent 
state in time. The derivation of time evolution is trivial if we restrict the evolution operator to 
be unitary. Non-trivial evolution occurs if the classically forbidden regions behind the horizon 
get exposed using a non-Hermitian term added to the Hamiltonian which facilitates entropy 
production [17,20]. The change in entropy of the black hole is carried away as Hawking flux 
which escapes to infinity. This time evolution is irreversible in nature and is very similar to such 
time evolutions in complex systems [34]. 

In the next section we briefly describe coherent states for gravity. We also review the concept 
of entanglement entropy in this section. The third section describes the origin of entropy using 
a density matrix. It includes a description of the tracing mechanism when the coherent state 
is defined using intertwiners at vertices. In the fourth section we describe time evolution and 
the relation of this entropy production process to physics of complex systems. The fifth section 
discusses the results and the open problems yet to be solved using this formalism. 

2 Coherent states in loop quantum gravity 

The coherent states are useful semiclassical states in a quantum theory. In case of the simple 
harmonic oscillator or in quantum electrodynamics, these can be formulated as the eigenstate 
of the annihilation operator. Explicitly 



where a is the annihilation operator for the positive energy modes of the theory and z is the 
eigenvalue. These are usually Gaussians, e.g. for the harmonic oscillator 

1 / (x-z)' 

tp = —j=^ exp 



V2^ V 2/i 

where z = xq ^ ipo is a point in the classical complexified phase space, h = /i/27r, and h is 
Planck's constant. Thus the coherent state is a Gaussian 'peaked' on the classical phase space 
point. In these coherent states the expectation values of operators are obtained as their classical 
values and trajectories also evolve along classical paths, preserving the coherence. 

The coherent states appear as the kernel of a transformation from the Hilbert space L^(M) 
to the Segal-Bergmann representation of the wave functions or H{C) n L^(C) [25]. Using the 
definition of the coherent state as a kernel Hall [25] generalised the coherent states to obtain 
those defined for a SU(2) Hilbert space. These appear as kernels in the transformation from 
the SU(2) Hilbert space to the intersection of the Hibert space defined in the complexified 
SU(2) phase space which incidentally is SL(2, C) space. These are therefore often named as 
'complexifier coherent states'. In a particular definition of the kernel, it is obtained as 

K{h,g) = e-^^6{h',h)h'^g, (2.1) 
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where h is a SU(2) element, and 5 is a SL(2,C) element, A is the Laplacian on the group 
manifold, t is a paraineter. The Hall coherent states can be directly applied to gravity as 
written in terms of LQG variables gravity is a SU(2) theory [42]. 

We begin by writing gravity in terms of LQG variables. The phase space of loop quantum 
gravity is obtained using a canonical formulation of gravity. The manifold is taken dif feomorphic 
to S xM. The S slices foliate the space-time, and using the ADM metric, the dynamical variables 
are identified as the intrinsic metric qab of the slice and the extrinsic curvature Kab with which 
the slices are embedded in the space-time. The intrinsic metric can be written in terms of the 
tangent space triads or soldering forms e„. The phase space for LQG is again a redefinition of 
these and described thus [21,22,41] 

Ai = Ti-l3Kabe''', Ef = ^idete)ef, (2.2) 

where e^ are the usual triads the 'square root' of the metric qab = e^e^, Kab is the extrinsic 
curvature, F^ are the associated spin connections to the triads, /? is the one parameter ambiguity, 
which is known as the Immirzi parameter. The quantisation of the Poisson algebra of these 
variables is done by smearing the connection along one dimensional edges e of length 6e of 
a graph T to get holonomies he{A). The triads are smeared in a set of 2-surface decomposition 
of the three dimensional spatial slice to get the corresponding momentum P/. (The momenta are 
also labeled by the edges as each edge intersects the corresponding 2 surface once by construction, 
or every edge has a corresponding unique two surface and vice versa.) The regularised LQG 
variables are 

he{A) = V e^p ( [ a] , pI=1\e^. (2.3) 



The algebra is then represented in a kinematic 'Hilbert space', in which the physical con- 
straints have been 'formally' realised. Once the phase space variables have been identified, one 
can write a coherent state for these [25], i.e. minimum uncertainty states peaked at classical 
values of he, Pj. for one edge of the graph. The tensor product of these coherent states for 
each edge of the graph gives the coherent state for the entire space-time. Using the Peter- Weyl 
theorem for the delta function 5{h, h') = X]i (^i + l)Xj {h'h~^^ , where Xj is the character of the 
jth representation of SU(2), the complexifier coherent state (2.1) can be written as 

4;\ge) = 5^(2i + l)e-*^-(^-+i)/2x,(5e/ie"')- (2.4) 

j 

These are coherent states peaked at the classical holonomy with a width controlled by the 
parameter t. A 'momentum representation' of this is obtained by defining a Fourier transform 
which is such that the 'momentum eigenstate' \jmn) gives 

{h\jmn) = -Kj{h)mn, Pe\jmn) = j{j + l)h\jmn). 

The \jmn) are thus the usual basis spin network states given by -v/2j~+T7rj(/ie)mn (normalised 
in the he basis), which is the mn element of the jth irreducible representation of the SU(2) mat- 



rix he, and Pe = a/P/P/. Taking the scalar product of the configuration coherent state (2.4) 

with the momentum eigenstate gives the momentum 'coherent state' as ipl = e~^^^^~^^'''^'Kj{ge)jmn, 
where 'Kj{ge)mn is the j'th irreducible representation of the SL(2,C) element ge- Thus the cohe- 
rent state in the momentum representation is 

W{9e)) = ^e-'^^^-^'y\{geUn\jmn). 

jmn 
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In the above ge is a complexified classical phase space element e*^ ^"^ ''^h'^, where P/'^^ and h'^ 
represent classical momenta and holonomy obtained by embedding the edge in the classical 
metric and T^ are SU(2) generators. The j is the quantum number of the SU(2) Casimir 
operator in that representation, and m, n represent azimuthal quantum numbers which run 
from —j, . . . ,j. Similarly, (2j + 1) x (2j + l)-dimensional representations of the 2x2 matrix g^ 
are denoted as 'Kj{ge)mn- The coherent state is precisely peaked with maximum probability at 
the /ig for the variable he as well as the classical momentum P^ for the variable P/. The 
fluctuations about the classical value are controlled by the parameter t (the semiclassicality 
parameter). This parameter is given by ti/a where Ip is Planck's constant and a is a dimensional 
constant which characterises the system. For the Schzwarzschild black hole a can be taken as 
proportional to the area of the horizon. The coherent state for an entire slice can be obtained 
by taking the tensor product of the coherent state for each edge which form a graph F, 

^r = nV'e- (2.5) 

e 

Thus we are considering a semiclassical state, which is a state such that expectation values 
of operators are closest to their classical values. The information of the classical phase space 
variables are encoded in the complexified SU(2) elements labeled as ge- The fluctuations over the 
classical values are controlled by the semiclassical parameter t. 

The density matrix which describes the entire black hole slice is obtained as 

where |^r) is the coherent state wavefunction for the entire slice, a tensor product of coherent 
state for each edge. These coherent states have been studied in various other forms [10,23]. In 
this paper we also discuss the SU(2) invariant coherent state described in [42] using intertwiners 
at the vertices. We will describe the detailed derivation of the intertwiners for a particular 
graph in the next section using this generalisation. We should mention that previously other 
semiclassical states had been considered in the framework LQG known as weave states [5], but 
none have been yet used to describe black hole entropy. Coherent states have been used to 
describe quantum cosmology and FLRW universes [32] . Relatively recent reviews on approaches 
to black hole entropy and corrections to black hole entropy in LQG can be found in [6,14,33]. 
A review of black hole entropy in LQG and recent description of two dimensional black hole 
evaporation can be found in [3] . Before we begin the discussion on derivation of entropy using 
density matrices, we briefly describe the concept of 'entanglement entropy' and its use in quan- 
tum field theory in curved space-time to describe entropy of black holes. Our derivation will be 
totally 'gravitational in origin' defined using the LQG phase space variables. 

3 Entanglement entropy 

In quantum field theory, we use the definition of entropy due to von Neumann. It is defined 
thus given a density matrix p 

S = — Tr{plnp). 

If the state is pure the entropy is zero, if the state is mixed the entropy is obtained as non-zero. 
This could be a system in which part of the wavefunction or the state of the system has been 
traced away. The tracing can be done for systems which have a product Hilbert spaces: e.g. 
H = Hi® Hq- The Hj is the Hilbert space which defines the 'internal space', Hq is the outside 
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Hilbert space. The basis states also thus appear in the product structure \i) (K> \j)- The states 
are said to be entangled if the coefficients of the basis states do not factorise, e.g. 

ij 

is entangled if 

dij / didj. (3.1) 

A density matrix defined thus 

P = I^>(V'I 

can have a partial trace performed on it, 

and the entropy computed using the above. We will now describe the 'entanglement entropy' of 
quantum fields in a black hole background. This shed light on why the entropy was proportional 
to the area of the horizon and not the volume of the black hole. A review on entanglement 
entropy of black holes can be found in [37]. Srednicki computed the entanglement entropy of 
a quantum field, where the field enclosed in a sphere of radius 'i?' is traced out. The entropy 
of the remaining space-time is proportional to the area of the sphere [38]. This was a very 
interesting calculation however doing a similar calculation for a black hole background did not 
give finite answers. The reason is that the 'metric' in the coordinates of an outside observer 
is singular at the horizon. The near horizon volume is infinite, and thus the entropy due to 
the quantum fields in this volume is also infinite, 't Hooft [40] introduced a mass dependent 
'cutoff as if a brickwall was steeling off the quantum fields, which satisfied a Dirichlet boundary 
condition at the wall. The presence of the brick wall could be used to obtain a finite value for 
the entropy. We shall briefly review the procedure here. 

3.1 The brickwall model 

The scalar field (f) is taken in the background of a Schwarzschild black hole whose metric in 
spherical coordinates is 

ds^ = -(l-'^\ dt" +(l- '^y\r^ + r^{d9^ + sm^9dcp'), (3.2) 

where M is the mass of the black hole, G = 1, c = 1. 

The scalar field's Klein-Gordon equation in this background is 



-9 



Using the inverse of the metric (3.2) one finds that the Klein-Gordon equation is 



- — j {-d^) <P + ^dr ir^ il - — \ drj (p + ^L^cj) - m^(P = 0, 

where L^ is the angular momentum operator. The above can be solved in the eikonal approxi- 
mation using 
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where 



k'= . ^' J fi--V'g^-A'('+i)->.»^ 



r 
r{r~2M) IV" f J ^ r^' 



As evident the wavenumber is quite big near the horizon, justifying the eikonal approximation 
(the phase dominates). Due to the Dirichlet boundary condition, (/> = at the wall, a particular 
value of the radius taken as r = 2M+h and an outer boundary L, the wavefunction's wavenumber 
gets quantised as integer multiples of vr 

nn = / kdr. 

J2M+h 

The total number of such eigenmodes with energy E is therefore the sum over angular momentum 
(each with 2/ + 1) degeneracy 

I {21 + l)dlTTn = g{E). 

Given this one can compute 



e-/3F 



n,l,l3 



where F is the free energy, and the righthand side is the thermal distribution, at a temperatu- 
re 1//3, such that 

9{E) 



ttPF = I dg{E) log (1 - e'^^) = - |" dE^^ _ ^ . 
Using the value of g{E) 



' r ,E f- drd-'-^V f (^' + ^).jE^-ri-^V'"+" 



-^ Jo J2M+h V r J J (e/3^-1) V V r J \ r^ 

(we use the approximation that the scalar fields are massless). F is then approximated in [40] 
as 

27r3 /2My 2 3 f°° dEE^ 



45h\i3j 97r Jq ^I^e _ i 
Using the definition of entropy 

S = J3{U-F), U = -t{^F) 

one obtains S = 4M^. This is indeed proportional to the area of the horizon, but one had to 
put the cutoff radius h = j20ttM ' ^ ^^^^ particular value to achieve the result. 

There are some observations about this derivation of entropy: (i) There is no entanglement 
of modes across the horizon of the QFT modes, what this represents is the entropy of a gas of 
scalar particles outside the horizon, (ii) Even though the cutoff is mass dependent the proper- 
distance of the brickwall from the horizon is a constant. The question still remains why it is 
this constant and not another one. 
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3.2 Scalar entanglement for a bifurcate horizon 

It is very interesting that this 'entropy of a gas of scalar particles' using the brick wall method 
could be placed in an 'entanglement entropy' context. In this formalism one traces over the 
'outside modes' to obtain a density matrix and the entropy of the resultant partially traced 
sector of the system. The details of this derivation can be found in [8] and references therein. 
The crux of the computation is hinged on the definition of the wavefunction of the black hole 
as a 'path-integral' in a space-time with a complex metric. The space-time is the complexified 
Kruskal extension of the Schwarzschild black hole. The metric in spherical coordinates for the 
Schwarzschild black hole is given in (3.2). One then defines the Kruskal coordinates as 

Tj _ _ -(t-r-2M ln(r/2M-l))/4M y _ (t+r+2Mln(r/2Af-l))/4M 

such that the Kruskal metric is 

ds2 = Z32^ exp-('-/2^) dUdV + r^dn\ 
r 

As evident the metric is extendible past r = 2M, the horizon, and the horizon is bifurcate, with 
V = coinciding with the past horizon and U = coinciding with the future horizon. There 
are two asymptotics, and the transformations U — )■ —U and V — t- —V map the two asymptotics. 
One could then define the Lorentzian metric as the real section of a complex metric, by defining 
the coordinates U and V w.r.t. a 'complex time' z = t + it. The r is the Euclidean time, and 
it is such that in the Euclidean section —2ttAI < t < +2ttM. The two asymptotics are then 
connected in some sense through this 'Euclidean' section. The complexified coordinates are 

TJ _ _ (iz-i2TTM+r+2Mlii{r/2M-l))/4M y _ {-iz+i2TxM+r+2M\ii{r /2M-i))/AM 

Clearly as the imaginary time r changes from — 27rM to 27rM, one travels from one asymptotics 
to another given by z'^ = 2ttM + it and z~ = —2ttM + it (— cxd < t < oo) in complex time. 

One defines a path-integral through the Euclidean throat which could be labeled as the 
propagator connecting the two asymptotics, with (j) = (p^ in one and (j) = (f)- in the other. 

Define 

*(</._,</>+) = (</._! exp(-W2)|</.+). (3.3) 

Thus the path-integral is over the Euclidean 'throat' connecting the two Lorentzian asymptotics. 
Given /3 = SvrM and if is a Hamiltonian which in the linearised approximation gives the usual 
Hamiltonian of matter fields including scalar fields. The computation in [8] is given for scalar 
Hamiltonian using heat kernel techniques. A density matrix is defined thus 



Plugging ^ from (3.3) and using j'P(/)_|_|</)_|_)(</)_|_| = 1 this reduces to 

p = exp(r)((/>_| exp(-87rM^)|0_), (3.4) 

where the prefactor is introduced to preserve JVcp^p = 1 or the trace condition. The entropy 
of this density matrix is then computed using S = — Tr(/>lnp). As the field is in one of the 
asymptotic regions, the Hamiltonian is obtained in the one-loop approximation as the due to the 
scalar fields solved in the classical black hole background. Due to the definition of the density 
matrix as the expectation value of the Hamiltonian at a temperature 1//3 (3.4), this calculation 
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can be eventually mapped to the computation of entropy of a gas of the scalar particles near 
the horizon and it is obtained as 

^~ 45 Lir-2Mr ^^-^^ 

Though this is proportional to the horizon area some aspects remain as the brickwall model. 
The integral in (3.5) is divergent and thus the entropy is infinite. The exact nature of the en- 
tanglement is not clear however in the way entangled states are defined in (3.1). 'Entanglement' 
entropy of scalar fields in flat space can be found where the entanglement is obvious as in (3.1) 
in the following references [11,38]. 

In the next section we will compute the 'entanglement' entropy of a black hole using 'non- 
perturbative' coherent states in a 'quantum gravity' formalism. The variables are 'regularised' 
gravitational degrees of freedom. The entropy computed is purely gravitational in origin. The 
nature of entanglement is also specified clearly exactly similar to the discussion of quantum 
mechanical entanglement in equation (3.1). The 'internal' Hilbert space is traced out and the 
entropy is computed using the density matrix defined in the 'outside' Hilbert space. The answer 
is finite and no cutoff is introduced. 

4 Gravitational entropy 

In the definition of the tensor product form of the coherent state (2.5), the total coherent 
state is a tensor product of coherent states at each edge. There is a SU(2) Hilbert space at 
each edge of the graph F. The total Hilbert space is thus H = ^H^. The question therefore 
naturally arises if the coherent state defined in the basis of these Hilbert spaces are entangled 
or simply independent? The structure, ^ = W^ ipe is such that the states appears independent 
in the tensor product Hilbert spaces, but can be factorised as per equation (3.1). However, 
the classical geometry is interwoven from edge to edge due to Einstein's equations. Thus even 
though not written in a manifest way; for two adjacent edges ei and 62, V'ei is related to ipe2 due 
to the classical solution at which these are constructed to be peaked. We label this entanglement 
due to classical geometry as classical correlations. We then compute the entanglement of edges 
outside the horizon with those inside the horizon. The coherent states as introduced in (2.5) are 
'gauge covariant'. The gauge SU(2) transformations act at the vertices and they are discussed in 
details in Appendix A. Introducing intertwiners at the vertices makes the coherent states 'gauge 
invariant'. They no longer transform due to the SU(2) transformations as the intertwiners map 
the vertices to the trivial representation of SU(2). These introduce further 'entanglement' of the 
coherent states of edges ending/beginning at the same vertex. These correlations arise as the 
type of spins assigned to the edges meeting at a vertex are restricted to ensure gauge invariance. 
These correlations introduced due to 'intertwiners ' are labeled as 'quantum correlations '. We thus 
identify two types of correlations, (i) The classical correlations which are due to the relation of 
the ge from one edge to the next already discussed in [18]. (ii) Quantum correlations due to the 
intertwiners which link the quantum spins of the edges meeting at a vertex. We discuss these 
correlations for the first time in this article. 

4.1 Classical correlations 

The classical correlations can be identified easily in a particular slicing of the classical geometry. 
We discuss a particular time slicing of the metric where the intrinsic curvature of the time slices 
is flat. One such metric which has the time slices as flat is the Lemaitre metric 

-|4/3 

rl/^{d9^ + sm^9d(l)'^). 



,9 ,9 dR 

ds^ = -dr^ + — - + 



l(n-r) 
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The rg = 2GM (in units of c = 1) and in the r = const slices one can define the induced metric 
in terms of a V coordinate defined as dr = dR/ [3/2rg(i? — Tc)] (r = Tc is a constant) on the 
slice. One gets the metric of the three slice to be 

dsl = dr^ + r^{de^ + sm^ed<tP'). (4.1) 

The entire curvature of the space-time metric is contained in the extrinsic curvature or K^y = 
i^drQixu tensor of the r = const slices. Now if there exists an apparent horizon somewhere in the 
above spatial slice, then that is located as a solution to the equation 

VaS'' + K^bS^S' -K = 0, 

where 5" {{a,b = 1,2,3) denote the spatial indices) is the normal to the horizon, Kab is the 
extrinsic curvature in the induced coordinates of the slice, and K is the trace of the extrinsic 
curvature. If the horizon is chosen to be the 2-sphere, then in the coordinates of (4.1), S"" = 
(1, 0, 0), the apparent horizon equation as a function of the metric reduces to 

Krr{l- q'') - K^<pq'^<^ - T^ - Keeq'' - Vl = 0. (4.2) 

In this article we provide a very simplified explanation of why this equation provides correlations 
across the horizon. Note that the first term of the equation disappears trivially as 1 = q^'^ for 
any point in the spatial slice and not just at the horizon in the classical metric. 

If we plug in the exact values of the above in terms of Schwarzschild variables, we find that the 
terms cancel groupwise at the horizon due to spherical symmetry. Thus the apparent horizon 
equation which respects spherical symmetry can be imposed simply as the equation in the 9 
variable or the (p variable 

This 'truncated' versions of the apparent horizon equation can be implemented on the co- 
herent states in the classical limit. Due to spherical symmetry, the other half of the horizon 
equation gets automatically satisfied. Note we are not using two different equations to replace 
the one equation in (4.2), but observing that due to spherical symmetry, ensuring that the equa- 
tion is satisfied for the 6 indexed variables ensures that the other half of the equation will also 
be zero. In other words the equations are proportional to each other. To regularise the equation 
we observe 

r^. = lq'\qe0,r) = -Iqeeiqf) = Iqee ^ {q'' (vi) - q'\v2)). 

In the above we wrote the Christoffel connection as a difference equation, vi and V2 are the 
vertices of the graph and 6e^ is the length of the radial edge connecting the two vertices. 

The above then translates to a equation for the truncated apparent horizon equation (4.3) as 

q'\v,^) = q'\v,^) [26,^Kie'M^) + l] , (4.4) 

where the fout vertices are outside the horizon and the Win are vertices within the horizon and 
Keeq^^ = K^e^j {K^ = Kabe^^). We then embed a graph with edges along the coordinate 
lines r, 0, (p of the spherical coordinate grid, and compute the regularised variables of (2.3). The 
details of this can be found in [18]. In terms of the regularised variables (4.4) appears as 



VPliVout)=Pl{v,^) 



2Tr (/.e-;(-in)T^/3^^^) PeAvin) + V 



(4.5) 
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where V is the volume operator and P^ = P^Pl- The volume operator appears due to the fact 
that it is the densitised operator y/qe'j which is regularised and y/q is regularised into the volume 
operator. The volume operator in spherically symmetric coordinates is written as PeuSe^ and 
thus it can be computed for this system. We also use a regularisation introduced in [18] for the iC^ 
operator which uses a 'derivative' w.r.t. the Immirzi parameter. Using the definition in (2.2), 
and the definition of holonomy one gets dhe/dp = (^fK^T^dx°')he- And thus multiplying 
with h~^T and taking trace identifies the appropriate component of the extrinsic curvature 
in the continuum limit. This regularisation works for the purposes of the computations of 
expectation values as verified in [20]. Note that constants which appear due to the specification 
of the graph edges (like 6e^) are suppressed in the above. And thus it is sufficient to take the 
following form of the apparent horizon equation 



VPl{Vout)=Pl{Vin) 



2Tr (/.e-;(^m)T^/3^^^) PeAvin) + V 



(4.6) 



We use the apparent horizon equation to solve for P^^ (t'out) in terms of Peg (um)- The apparent 
horizon equation is further supplemented with a restriction that the area induced on the horizon 
using the radial edges sums up to the total area of the horizon. Which translated into the phase 
space operators means 



E^- 



Ah. (4.7) 



ej/ 



These (4.6), (4.7) therefore specify the horizon information. We then use these to perform 
the tracing mechanism on the density matrix and compute entropy. Instead of computing the 
entropy with the entire density matrix, we concentrate on a local set of vertices surrounding 
the horizon. As we saw, the apparent horizon correlates 'angular edges' immediately inside and 
outside the horizon. We label this element pi°'=^' and isolate it from the entire density matrix. 

total outside local inside 

P " P P ) 

where p^"'^^^ covers a band of vertices surrounding the horizon one set on a sphere at radius 
^g — ^sr/"^ and one set on a sphere at radius rg + 6e^/2 within the horizon, as described in [18], 
and in the Fig. 1 enclosed. This local density matrix and the correlations due to the apparent 
horizon equation (4.6) was used to derive entropy [18]. 

We then further concentrate on two vertices Uout outside the horizon, vin inside the hori- 
zon which share a radial edge en and write a density matrix for that p'°'='^^. At each vertex 
■Vout,in there are two angular edges corresponding to the 9, coordinate lines which are in- 
going and two angular edges which are outgoing. These angular edges, 4 in number at each 
vertex have their classical ge correlated with those meeting at the other vertex. Labeling 

4 
Yl Yl'4'jo^(9e{Vout)[ge{v■m)])mo,no,i^jH{9eH)m„nH 

{jo^}{jH}{jIi},{j'oi}{j'H}ijI^}^' = ^ 
4 4 

X '[l'^JHi9e{vin))mnnn JJ \{JOi}{JH}{jli}){{Joi}{j'H}{Jli}\ 



local 



4 4 

X 

4 = 1 j = l 



]Xi'j'^^{9e{Vont)[9e{Vm)])ni'Qy^^i^j'jj{9eH)m'^n'g]\4^j'j^{9e{Vm))m'j^n^^ 



JOiT^Oi'^Oi label the spins of the edges at the outside vertex, jumnnii label the spins at the 
inner vertex and the jh, mn, nn label the spins on the radial edge connecting the two vertices. 
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Horizon 
Figure 1. Graph at the horizon. 

We us {j} to label all the three indices jmn. The tracing over the modes inside the horizon 
gives in the t — t- limit 



local 
rreduced 



^ Yl'4'jo^(9e{Vout){9e{Vm)))mo,no,'iPjHi9eH)mHnH 

{J0r}{JH}{jli}{J0i}{fH} * = 1 
4 

xY[^{Pe,{Vin)Jnt)\{jo^}{jH}){{j'o^}{j'H}\i^fJ9e,) 



' n' 



i=l 

4 

X n ^j'oS9eiVout)[geiVin)])m'o,n'oi- 



In the t — )• limit, the diagonal terms dominate. The density matrix collapses to a set of 
non-zero elements [42]. This can be directly seen from 



-*^'(^'+^Vj(5-)„„7rj(fif), 



/ /I /\^ '' i\y Jmn'' 'i\y jmn 

^3/2 



4:y/Trp 



g-ip|m/i-(^P)3/p|g-jp|"/j-(^p)3/p|g-[(j+l/2)-p]2/i_ 



{U 



In the t — )• limit, this gets to be peaked like a delta function at the values P = {j + l/2)t, 
mt,nt = ^P^j^P^ are the left and right invariant momentum operators. As in [18] the joi gets 
fixed to be a certain number due to the apparent horizon equation (4.6), and the correspon- 
ding moi, noi were also taken as fixed numbers using a 'gauge fixed' version of the apparent 
horizon equation. However this is slightly arbitrary as the classical data is gauge invariant. We 
expect that when all the constraints including the Gauss constraints are solved this will emerge 
as obvious. 

The individual components of ^^'P^q ^ fixes the moiriQi to unique numbers. Thus the outer 
classical data Jq^, rrfov ''^Oi ^^^-^e to be specified completely to specify the 'outer state'. In other 
words we can break the degeneracy of the density matrix with different moinoi by making the 
operator measurements. 

Now we come to a subtle point I had not discussed before, and that is the fact that the part 
of the radial edge en (see Fig. 1) is hidden behind the horizon. It thus makes sense to trace over 
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that. To implement that, we therefore divide the horizon edge en into two further edges eni 
and eH2- We build the basis state for this using the fact that 

\jmn) = i/2j + llTj{he)mn\0) 

= ^ \/2j + lTTj{he^)mk'^j{he2)kn\0) = ^ , \jmk)\jkn), 

k fc V^J + i 

where |0) is the vacuum state. Splitting the horizon edge, and tracing over the internal half 
gives 



Tr 62] jmn){j'mn \ = > — -, \j'mk){j'k'n\jkn){j'm'k'\ 

k'kn 



= Y^ ^-T^|jmA;)(jm'/i;|(^„„/. (4.9) 

k 

Similarly, for the irreducible representation components of the 'wavefunction' 

^i \9e)inn'^j \9e)m'nj 

which multiplies the basis states in the density matrix, one can perform such a splitting of the 
holonomy as 

„ -JT'P'f^ -JT'P\ ^ _ ^ 
i/e — '-■ "'e — c; /(,ei"'e2 — i/ei "-62 ) 

and thus 

/ , ^j (5ei )mk'^j{he2)knT^j {Qci )m'k'T^j{he2)k'n- 
kk' 

The sum over n then gives 

■^j{gei)mkTi-jige)m'k = '^ {f^''^' ^' ) mm' ' (^-^O) 

This and then its complex conjugate gives the delta function in the t — t- limit peaked at Pjj- 
The contribution from the horizon edge to the density matrix is 



— Y X] ^iJHt, Pe„)\JH'mHkH){JHmHkH\. 



mnkH 

This explains the origin of the 1/{2jh + 1) factor anticipated in [18]. 

Thus the reduced density matrix 'describing the vertex outside the horizon' is 

mnriH 

4 

X n l^eo»/P30i ^P30i)\PeH'mHnH){PeH^HnH\{Peo^,^P30i ^P30i\- 
i=l 

And the tensor product density matrix becomes 

1 

n l-f'eo./psOi ^PZOi)\PeH^HnH){PeH'mHnH\{Peo^^^V?,Oi ^P30i 



local = "TT V ^ 

^reduced H /_^ 2jfi + l 

■"out muinH 



4 
X 

=1 
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Figure 2. Spins at the vertices. 



the above density matrix is clearly diagonal. If we then compute 



T^^l' Jocal 1 local \ 
-'-'^ \Preduced rreduced/ 



4n 



+ corrections, 



where /3 is a prefactor [19]. The exact computation is not obvious, but if we assume that 
every jh is the same, then this gives - Tr (/3[°duced l'^ /'reduced) = E^^.t ln(2j// + l) = Ah ln{2JH + 
l)/[ll{JH + 1/2)] where we used the constraint (4.7) such that the Y.{JH + 1/2) = Ah/II- For 
different distribution of spins one has to use combinatorics [18, 19]. 

Thus to leading order the entropy as computed from the density matrix is indeed Bekenstein- 
Hawking modulo the (5 which is a function of the Immirzi parameter and other details of the 
combinatorics of the graph can be set to 1. 



4.2 Quantum correlations 

We now discuss the correlations across the horizon due to intertwiners placed at the ver- 
tices. As evident the coherent states are defined in a SU(2) Hilbert space. For the holonomy 
variables, the SU(2) gauge transformations act at the beginning and end point of the edges. 
hf. — )• g~^{vi)he{vi,V2)g{v2) ■ These changes will change the spin network state. To make the 
state gauge invariant one imposes intertwiners at the vertices. These are structures which map 
the spin networks at the vertex to the trivial representation. Thus the state does not transform 
due to the intertwiners. We concentrate on two vertices (suppress one dimension) and take each 
vertex to be four valent at the horizon and observe the role of the intertwiners. From the Fig. 2 
two of the edges are ingoing, the other two are outgoing at each vertex. The gauge transforma- 
tions act on the holonomies and hence the basis states on each edge at the vertices. Integrating 
(group averaging) over the gauge group at the vertices, one ensures gauge invariance [21,22,41]. 
This solves the Gauss constraint. We find that when done carefully this gives us correlations 
between the coherent state at inner and outer vertices which are carried by the horizon edge 
which connects the two vertices. This is therefore a 'quantum correlation'. We describe some 
of the details of the calculations here. Let us take the two vertices as win and Vout and the 
spin labels j/i, jn-, J13, Jh at the inner vertex and joi, J02, io3; JH at the outer vertex. Let 
us assume that at the inner vertex jn, jh are outgoing at that vertex, or they 'begin' at vin 
and the J12, J13 are ingoing at the vertex or have their end point at that vertex. Similarly at 
the outer vertex, jh and jo3 are ingoing or end at fout and joi, J02 are outgoing and begin 
at the vertex t'out- To introduce the intertwiners we observe the basis spin network functions 
and their transformations. The basis spin network states \jmn) at the inner vertex in the he 
representation are given as {\/2j + 1 factor is not shown but it exists in each basis state for 
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normalisation of the inner product) 

T^jjl l"'jm/inji7^j/2 ("'/m72"/2'''"j/3 i'^jmi^nj^'^JH \'^)mHnH , 

the gauge action U{g) at the vertex is such that 



'mi2ni2 -'JJ V -^ ^rnisnis Jti \^ jmnriH- 



Upon group averaging one obtains a projector onto the gauge invariant Hilbert space which 
written using the 3jm symbols is thus (see Appendix A for further details) 



^-^ \mn ruH rukj \pn pn rikj 



{jk}piiqi2qi3PH 



^ , J 12 J 13 Jk \ I J 12 J 13 Jk 
X 



v"/2 n/3 muj \qn qi3 nk^ 

We can then collect the basis vectors and the 3jm symbols contracted with them as a |in) basis 

^v^/^i/i JH 3k\fin ji3 jk 
j-^ \mn ruH mt) Kjin nj^ nik^ 

If we take 'ipj{ge)mn = (m^\ g~^ '' '^j{9e)mn, the coherent state wavefunction at the inner 
vertex is 

Z^ii'^luKyejJmum^y^^^ ^^ ^^^J ^^n n/3 mfcy/ ' ' 

{3} «=1 

Note that this way, only the indices which begin or end at the fin are gauge invariant, the 
other end of the edges remains 'free' and can be attached to another vertex in an invariant way. 
Exactly in the same way one can obtain the gauge invariant wavefunction at the outer vertex, 
comprising of the 3jm symbols and the corresponding wavefunctions. And therefore one can 
define the intertwiners at the outer vertex but notice immediately the horizon edge is shared 
by both the outer vertex and the innner vertex, and thus, the correct way to write the coherent 
state for the edges meeting at the two vertices is this 

|^(^'out,Win)> =Yj^i{J0i},{3H},{3ii})\ont)\H)\m), 
{J} 

where \H) is the horizon state for the horizon edge crossing the horizon (see Fig. 1). Clearly the 
coefficient of the basis states do not factorise and can be identified as 

o ~ 

i=l 'At ^ ^'^ 

Jk 

( jk' J03 Jh\ I / X 

Xnik' no3 nkj 



^ Y.W'i'Ji^^S. 



J II JH Jk \ I J 12 J 13 Jk 

eiiJmjinji 



. . ^mn niH rukJ \nn nis nik 
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The density matrix is then 

Pinv = |V'Kut,'Vin))(V'(^out,'Uin)|- 

This in a very sure way establishes that correlations exist across the horizon surface, the 
question however is that is Bekenstein-Hawking entropy recovered from the tracing mechanism 
as previously? Ab initio it is not expected that this will happen, as these gauge invariant 
coherent states do not have the appropriate peakedness properties of the gauge non-invariant 
coherent states [42]. We simply implement the tracing of the inner edges to find the reduced 
density matrix at this stage and demonstrate the entanglement. 

The reduced density matrix is easily obtained, using the trace over the inner basis states and 
some 3jm symbol identities (see Appendix A). The matrix elements turn out to be 

Yl * (fe}' ij^}' fe})'^({j/J, {J'h}, {j0^}) 
{jli} 

,T,a ■ \ f ■ W^J.ff ■' \ f ■' W \^ (2jfc + l)S{jliJHJk)S{jkjl2Jl3) 

= ^(uoi|,ml)*({jHl,{JoJ) 2^ ^2j + 1) 

{ju},jk 
3 3 

i=l i=l 



311 JH 3k \ / 3k 312 313 \ / 311 3H 3k \ 3k J/2 J/3 



mil rnu rukj \mk n/2 rij^j \m'j^ m'^j m'/^J ^m'^ n'j2 n'j^ 



(4.12) 



Whereas it is not obvious what this will yield for the entropy which is 5bh = ~Tr(/9ln/9), 
one can use the fact that again in the limit that t — )■ 0, the product of wavefunctions will assume 
the form in equation (4.8) 

Y^,T, /■ • • A,T,^- •/ •/ ^ ,T,/- ■ \if,fi -I yS:r(.yk + 'i^)S{jiiJHJk)S{jkjl2Jl3) 
Z^"^ {30i,3H,3li)'^{3li,3H,3oi) = ^{30i, 3h)^{3h^ 3oi) Z^ ^^^ +1) 

{ji} ik 

jfi 3H 3k\ f 3k 3% J/3^ (fk JH Jk\ f 3k J/2 Jl3 
mj\ TiiH rukJ \mk nf^ "'/s/ V^fi "^// '^'kJ V^'k "-/i "-/s 

The above reduced density matrix is precisely of the form expected to yield the entropy, but 
the 3jm symbols are non-trivial and one can use some asymptotic form given in [35] for large j 
to estimate the form of the density matrix. Details of this calculations and further computations 
of the entropy will appear in [15]. This is the absolute proof of the fact that the coherent state 
introduced to describe the black hole is entangled across the horizon and does not factorise into 
inside and outside coefficients. There has been work in observing entanglement in semiclassical 
states independently as in [30]. 

5 Time evolution 

In this section we discuss time evolution of the system. Fortunately as the system is semiclassical, 
a natural choice of time is the asymtotic time coordinate. Using the Lemaitre slice and evolving 
the state in Lemaitre time would not serve the purpose as that is the frame of an infalling 
observer. We use the frame of an observer stationed outside the horizon to evolve the slice in 
her time as at the asymptotics this coincides with Minkowski time. Again we concentrate on 
the evolution of a local patch surrounding the horizon. Normally the vertices within the horizon 
are classically forbidden to the outside observer. However we allow for access to those vertices 
in the 'quantum mechanical' evolution. The 'Hamiltonian' which is used to evolve the horizon 
is the Brown and York quasi-local energy. 
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5.1 Semi classical Hamiltonian 

Quasi-local Hamiltonian is defined as the energy enclosed in a finite space. This is obtained using 
a particular definition [13], as the 'surface' integral of the extrinsic curvature with which the 
surface is embedded in three space. This generates time translation in the timelines orthogonal 
to the spatial slice. In our case, we take the bounding surface to be the horizon and the quasilocal 
energy is given by the surface term [13, 28] 

H = — d Xy/ak, 



where k is the extrinsic curvature with which the 2-surface, which in this case is the horizon S^ 
embedded in the spatial 3-slices, and a is the determinant of the two metric a^^^ defined on the 
2-surface. This 'quasilocal energy' is measured with reference to a background metric. Thus 
H = H — Hq where Hq is the energy in the background. We concentrate on the physics observed 
by an observer stationed at a r = const sphere. 
The metric in static r = const observer's frame is 

ds^ = -fdt^ + r^ {de^ + sin^ Odcp^) , 



/ = ^JY^^V^fr^ where Vg is the Schwarzschild radius. If we take n^ to be the space-like vector, 
normal to the 2-surface, then the extrinsic curvature is given by k^^y = o"°Vc^?^t/ and the trace is 
obviously 

A; = V"n„. (5.1) 

In the special slicing of the stationary observer the normal to the horizon 2-surface is given by 
(0, /(r), 0, 0). However, we built the density matrix on the Lemaitre slice. The Lemaitre and the 
Schwarzschild observer's coordinates are related by the following coordinate transformations 

"" dr = {dR-dT), dt= -^-{dr- f'dR) , /' = ^. 

The r = const cylinder of the Schwarzschild coordinate corresponds to dR = dr of the 
Lemaitre coordinates, and for these dt = dr. Thus unit translation in the t coordinate coincides 
with unit translation in the r coordinate. Further, the intersection of the r = const cylinder with 
a t = const surface coincides with the intersection of r = const and the r = const surface. Thus 
in the initial slice, the QLE Hamiltonian can be written in terms of the Christoffel symbols of 
the Lemaitre slice. Even though the generic transformation of the Christoffel symbols from one 
coordinate frame to the other is inhomogeneous, in this example, the transformation is trivial. 
The trace of the extrinsic curvature (5.1) then written using Christoffel symbols turns out to be 



2k 



-^'-^-9 



dr dr 



fir) -Ho. (5.2) 



The reference frames' quasilocal energy Hq is a number, it just defines the zero point Hamilto- 
nian. Thus, we replace the classical expressions by operators evaluated at the r = const slice. 
In the first approximation we simply take the /(r) as classical y^l — rg/r = ^J^^J^h^g = e, as 
this arises due to the coordinate transformation and the norm of the vector rir in the previous 
frame. In the re-writing of (5.2) in regularised LQG variables the Hamiltonian appears rather 
complicated. 

One can rewrite these in a much simpler form, using the apparent horizon equation. Since 
the Hamiltonian is an integral over the horizon, the variables will satisfy the apparent horizon 
equation (4.2) up to quantum fluctuations. Thus the Hamiltonian operator is then re-written as 

i^horizon = - / d9d(P,/Mm[^ee^'^ + Kj^e^"^] , 

K .1 
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where we have used the classical apparent horizon equation (4.2) (with q"' = 1). The regularised 
LQG version of this is 



_ Cae sr-^ 

-^horizon — n r / ^ ^'^ 

ZKOegSeg 

t^out 



K'T^^-^K 



Pl+h.c. + ie^<P), (5.3) 



where C consists of some dimensionless constants Seg is the 2-dimensional area bit over which E^ 
is smeared, a is a dimensionfull constant which appears to get the P/ dimension less. 6eg is 
the length for the angular edge eg over which the gauge connection is integrated to obtain the 
holonomy. The sum over Uout is the set of vertices immediately outside the horizon. The (5.3) 
can then be lifted to an operator. However, this Hermitian operator gives no change to horizon 
entropy as the coherent state is evolved from one time slice to the next. Surprisingly a non- 
Hermitian term evolves the slices non-trivially. These anti-Hermitian terms can be easily found 
in this formalism if we allow for the vertices within the horizon to be exposed to the outside 
observer. In this case as the region is within the classical horizon, the norm of the Killing vector 
is negative, and n^ has components which are imaginary. The e — )• ib/,e. One sees why the 
inner vertices are relevant in a semiclassical evolution even if they are classically forbidden: In 
quantum mechanics tunneling is a well known phenomena, and inclusion of the inner vertices 
facilitates the process. 

Thus the contribution of 'inner vertices' Vjn is added to the Hamiltonian and it gives the 
anti-Hermitian contribution to the Hamiltonian. Thus -ff horizon = M^^,, H^^. + V,,. H^.]. 

di L * — ^ t^out (JUL i. — '(/in 111 J 

The regularised Hamiltonian is not Hermitian, and the evolution equation is 

ih^ = Hp-pH\ (5.4) 

Thus one takes the evolved slice to have a density matrix 

The entropy computed using this density matrix is different from the entropy in the previous 
slice and the difference of entropy can be seen as obtained using a simplified model. 

Let the density matrix be defined for a system whose states are given in the tensor product 
Hilbert spaces Hi ® H2 and given by 

\ip) = ^dij\i)\j), 

where \i) is the basis in Hi and \j) is the basis in H2 and dtj are the non-factorisable coefficients 
of the wavefunction in this basis. Let us label the wavefunction at time i = to be given by the 
coefficients d^-. The density matrix is 

iji'j' 

The reduced density matrix if one traces over H2 is 

ii' j 

We now evolve the system using a Hamiltonian which has the matrix elements Hijirj'\i)\j){j'\{i'\, 
we assume that the Hamiltonian does not factorise, that is there exists interaction terms between 
the two Hilbert spaces. The evolution equation is 

if^Q- = [H,p], 
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which in this particular basis gives the density matrix elements at an infinitesimally nearby slice 
to be 



^i'j'^ij ~ d'i'j'^ij ~ J.^'^ 



/ , {Hijkidj,idi*j, - dijd,^*iHkii'j') 



kl 



Thus we evolve the 'unreduced' density matrix and then trace over the H2 in the evolved slice. 
The reduced density matrix in the evolved slice is 



E^S*^& = E^^4-> 



2Z {^i3kidlid^*j - d°jdl*iHkii'j] 
kij 



This gives 



pSr ^pO_ ^^^^^ 

where A represents the commutator. Clearly the entropy in the evolved slice evaluated as 
6*3^ = — Tr(/9ln/?) can be found as 

4h = SU+'-5t[TvA In pO + Tr pV° "^A] . 

Given the definition of An, one gets 

A-ii = 2^ [pijkiHkiij - HijkiPkiij] ■ 
jki 

In case both the Hamiltonian and the density operator are Hermitian, one obtains 
^Ajj = 2iImTr(p°i?). 



This is clearly calculable, and gives the change in entropy AS'bh- The Inp^ term yields 
corrections, and we ignore it in the first approximation. However if the Hamiltonian is non- 
Hermitian one gets 



A5BH = ^Tr[i7/-pOFt] 



(5.5) 



instead of the commutator. 

In case of the gravitational Hamiltonian, the above can be computed using a U(l) projection 
thus. Let us take the U(l) case to make the calculations easier and observe the action of the 
Hamiltonian on the evolution of the coherent state. The spin network states are replaced by 
|n) = e'"^, < (^ < 27r, n is an integer and the coherent states are 



i;\g,) = Y,e- 



(tn2)/2 m(xe-Jpe)„-mC 



due = e^"'^^" ^P"' is the complexified phase space element in the 'nth' representation. 
The quasi local energy operator (5.3) also takes the simplified form 



H, 



U(l) 



horizon 



1 - 5 - 1 dh^^ - 

--C'Lh~^l3 — hePe + -C'LPel3-^he. 



(5.6) 
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The prefactors have been clubbed into C . The U(l) coherent states are eigenstates of an 
annihilation operator defined thus: 

The holonomy operator can thus be written as 



e-"'e-f^g. 



And the derivative w.r.t. Immirzi parameter of the holonomy which appears in the definition of 
the Hamiltonian replaced by 



/3 



dhe 



-t/2 



4-'^--"-^i 



-t/2 



p,e-P^g, + e-Pf^^ 



The dependence of the operator p on the Immirzi parameter is known (2.2), and thus we 
could evaluate the derivative {I3dppe{li) = (3dp{pe{l)/(3) = —pe{/3)). 
The term 



0[rU(l) 



^^ \P -"horizon^ 



(5.7) 



is then computable. Let us take the first term of (5.6) and find (5.7). As p^ = \ip){ip\, (5.7) 
gives simply (we drop the 'e' label for brevity) 



1^, 



u(i) I 

horizon I 



lc'ih-'l3-^hm + {^P\h.c. 



We then concentrate on the 2nd term of the above 



iC'{^l^\g^e 


-t/2g-j5g-t/2 


P e-^g + e-^13^ 


pli') 


+ (^|h.c.| 


iC'e-'g* 


{iple-Ppe-Pp 


^)5 + (V|e-2^/3|^ 


p\i') 


+ (V'lh.c. 



{^\e-'^/3^p[ 



'^'^% ldH9'W)mm) 



M9')(3j^{i^\e-'^m{^' 



where we have used the fact that coherent states resolve unity. It can be shown that the 
expectation value of the operators in the t — )• collapses the integral to g' = g point [42]. Thus 
one obtains from the above 



Tr{p'H^^J=-lc'e-^/' 



p + g*e-^Pp 



d(3 



p + h.c. = C'l3-7^p, 



(5.8) 



which is real, and thus 
A5bh = 



(this is actually the classical quasi local energy as it should be from Tr(p°//horizon))- 

This is obvious, as the way the Hamiltonian is defined, this is simply a function of the 
Hilbert space outside the horizon, and the matrix elements of this will not yield anything new. 
We approximated the horizon sphere by summing over vi vertices immediately outside the 
horizon. We could do the same by summing over V2 vertices immediately within the horizon. 
For the Lemaitre slice, the metric is smooth at the horizon, and one can take the 'quantum 
operators' evaluated at the vertex V2- In this case however, as the region is within the classical 
horizon, the norm of the Killing vector is negative, and n^ has components which are imaginary. 
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The e — )■ itte. Thus /:^horizon = ^[Si,! ^vi + X^^j -f^i'2]- ^^ ^^^ evaluation of the quasi local 
energy term, the energy would emerge correct in the Se^ — ?• limit as e — )• 0. The regularised 
Hamiltonian is not Hermitian, and the evolution equation is (5.4). 

The same type of expectation or trace is obtained of the imaginary term (or V2) is obtained 
as in (5.8), but this does not cancel in (5.5), and 

The 'rate of change' of entropy is thus 



A5bh = Ty2/5^>77 P 



if dp 



(we extracted the k from C" to get ti and rewrote the rest of the constants as C) . 

Thus there is a net change in entropy, but, to see if this is Hawking radiation, we have to 
couple matter to the system. 

In the case of the black hole system, due to spherical symmetry, in the classical limit, from 
the regularisation of the extrinsic curvature components Kq , which appear in the Brown and 
York quasi-local energy and its LQG regularisation (5.3), the change of entropy is 



A5bh = T-^J;/3 



P 
y 



d/3 '"> ^ 8(3 '^ 



1 


f ^ 


1 


— 


' *Ee = 




a 




a 



Due to the nature of the classical metric, the P/ can be gauge fixed to have one non-zero 
component in the internal directions [20]. Let that be some fixed index r, then P^g . = 5JP/ . 
In the computation of the Hamiltonian then it can be shown that if the holonomy is assumed 
to be he = e*-^ then, the r component of x^ contributes to the Hamiltonian, and this is what 
we have set as Xeg ^ in the above. From [16] one finds the Xee ^^^ the derivative of that at the 
horizon is a function of \/rg/r which gives a constant. The Peg are the regularised densitised 
triads at the horizon and thus they give area of the two surfaces comprising the decomposition 
of the three sphere in the angular dimensions. As it is evident, 

*r sin^. 

Thus they will be proportional to rg sin 9 at the horizon. Similar analysis is done for the (p sector 

, 2C'e6T V- , . o 
ASbh = ± — j2 — 2^ dAyfirg, 

where dA^ is the area of the element of the two surface at the vertex sm.9 and C" has been 
redefined to accomodate the regularisation constants which appear in Peg- 

Thus using the fact that Th = l/invg one could define the 'rate of entropy change' ASbr/St 
as 

ASbh oc ±--. 

Note that this is change in entropy due to emission of one Hawking quanta of a particular 
frequency. It should not be confused with entire black hole decay rate. 
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5.2 Hawking radiation 

The fact that entropy change occurs due to Hamiltonian evolution which includes a non- 
Hermitian term is very interesting. Hawking radiation however is a flux of particles emerging 
from the horizon, and thus we have to trace the origin of such a flux. In the previous section we 
found that as the system evolves in time, the horizon fluctuates and the area decreases. But is 
this Hawking radiation? Adding matter to a 'coherent state' description of semiclassical gravity 
has been discussed [36]. Thus, given a massless scalar field Lagrangian coupled to gravity, whose 
Hamiltonian is given by 

^2 



Hsc = I d^x 



the 'gravity' in the Hamiltonian can be regularised in terms of the he, Pj. operators in the 
coherent state formalism. The integral over the three volume gets converted to a sum over the 
vertices dotting the region. Thus 

V 

This Hamiltonian is an operator, and one evaluates an expectation value of the Hamiltonian 
in the density matrix by computing the trace of the product of the density matrix and the 
Hamiltonian 

Tr(p^//i). 

This Hamiltonian and the density matrix are then both evolved according to the time-like 
observers frame. One gets 

Thus one can compute the Hamiltonian in a slice infinitesimally close to the previous slice t-\-5t. 
Using that, 

Tr (p-+^-F;+^-) - TV {p-H-^) = -{6TfTr{[H,p-][H,H-^]}. 

It is very clear thus that the order 6t terms are zero for this. However, allowing for the 
non-unitary evolution using the non-Hermitian Hamiltonian, the 6t terms survive. In fact the 
terms are 

Tr {p^+'^H:+'-) - Tr [p^Hl) 

= -^ TV [{Hp - pH^)H,,] - ^ Tr [p{HH,, - H,,H^)] . (5.9) 

In the exact Hilbert space of the matter gravity system, the states are not exactly in the tensor 
product form, however, in the first approximation we take the states of the matter gravity 
system to be Y\^ 1^(5)) (X" \uj) where |a;) is a matter state with particles of energy uj. At this step 
for simplicity we assume that the outer and inner vertex scalar state is at the same energy uj. It 
is obvious that the first term in (5.9) gives the change in entropy as per (5.5) times the matter 
Hamiltonian's eigenvalue ASBHi^- The 'rate' of particle creation thus has the form 
euj 

where Th is the Hawking temperature for the signs +(— )te and negative (positive) uj. The e 
can be interpreted as the redshift factor. If a n particle matter state is taken, then the number 
will change accordingly. Thus a flux of particles indeed emerges from the horizon. To obtain 
a complete finite flux one has to compute the flux production over a finite number of such time 
steps. One can see that eventually a exponential will emerge which is the required Boltzmann 
factor. Further details can be found in [20]. 
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5.3 Irreversible systems 

It has been known from the time of Boltzmann that entropic systems usually have irreversible 
time evolution. The Boltzmann H-theorem uses a particular collision term in the evolution 
equation whose behavior gives rise to this irreversible flow. In complex systems physicists have 
been using different techniques to create irreversibility in the dynamics. In [34] this flow of 
a complex system is formulated as 

where the usual Liouville operator has been written as broken into two pieces, one ^^ which is 
the reversible part and the other <&^ which is irreversible and non-Hermitian. The ^^ creates the 
entropy production during the flow of the system. Thus it is remarkable that in the derivation 
of the time flow in presence of the horizon, precisely such a splitting occurs in the time-evolution 
equation. It seems quantum gravity should be formulated such that its dynamics should not 
be restricted to unitary operators. A complete re-formulation of quantum gravity using the 
language of complex systems and verification of some of the existing results is work in progress. 
This will facilitate the description of the system to describe black hole evaporation. Notice that 
as the non-Hermitian component of the Hamiltonian contributes from within the horizon, if all 
the mass of the black hole is radiated away, the term within the horizon will also disappear, and 
the evolution will become unitary again. However this as of yet speculative, a stable Planck size 
remnant might also form in the process. 

6 Conclusion 

Thus in this article we discussed coherent states for non-rotating black holes. We then showed 
that these states can be used to explain the origin of entropy due to classical correlations and then 
gave a brief introduction to a computation of 'quantum correlations' using SU(2) intertwiners at 
the vertices. I then derived a time evolution with a non-Hermitian Hamiltonian which generated 
entropy and this seems to create a flux of matter which escapes the horizon. This semiclassical 
derivation is similar to the process of entropy production observed in complex systems. However, 
the non-unitary nature of the evolution equation might be due to the use of semiclassical 'time' 
and an evolution operator which generates evolution in that time. One has to further work 
with a 'quantum Hamiltonian' to verify these results. Of course this brings us to the problem 
of identifying time in quantum gravity, and the evolution in physical time remains a project 
of research. There is promise from the reduced 'spherically symmetric' solution for the true 
Hamiltonian which can be found in the following papers [1, 12,24]. Other examples of reduced 
phase space spherically symmetric quantisation and studies of Hawking radiation can be found 
in [29,31]. For the time being, reformulating quantum gravity in an attempt to answer all the 
questions is work in progress. 

A Intertwiners 

To build the 'gauge invariant coherent states' we use the techniques of group averaging, the 
details of which can be found in reviews of loop quantum gravity [21,22,41]. For this the graph 
and its vertices are specified and the intertwiners at each vertex introduce the gauge invariance. 
For us, it is sufficient to take a graph which is comprised of vertices immediately outside the 
horizon and those immediately inside the horizon, connected by a radial edge. If we suppress one 
angular dimension of the spherical horizon, and open up this graph, it will appear as a 'ladder' 
structure. 
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If we isolate one rung of the graph and concentrate on that one obtains these two vertices f out 
connected to Vm with four edges meeting at each vertex (see Fig. 2). From the figure it is clear 
that two of the edges are ingoing and two are outgoing at each vertex. To obtain the intertwiners 
at each of the vertex, we write the basis states. At vertex Win the basis states are (suppressing 
the \/2j + 1 factor which makes the inner product 1) 

'^jli \''')mjinii'^jr2 \'^)mi2ni2''^ji:i K'^jmisniz'^JH K'^jmHriH 

Implementing a Gauge transformation at fjn one obtains 

~ ^i/i y9)miipji'^jji Wpiinji'^ji2 K'^)mj2qi2'^jl2 \9 ) qj2nj2 
^ '^jlsW^lsQls'^JiaKd )qj^nis''^iH\9)mHPH'^JHWpHnH- 

Integrating over the group element using the measure J'D{g) then gives the 3jm symbols. This 
then defines the projector onto the gauge invariant space. If we isolate the integrals this is 



j Vg TTj^^{g)mj,pj,7rj^2{9 ^)g,,„,,^i/3(5 ^)g^^^nr,''^„i9) 
Inserting a delta function using Peter-Weyl theorem 

^ (2jfc + l)7r,,(/i)„,„,7r,,(/.'-i)^^„^ = 6{hh"'] 

One obtains 



m„pH- 



jknikUk 

X J ^9T^h (5 "')^,n,^j,2 (9'"'),^^„^^vr,i3(5'-i)g,3„,3, (A.l) 

these integrate to 3jm coefficients as per the definition [43] 

/"'n„_ I „\ _ i„\ _ („\ _ ( 2\ h h\ ( h h h 



^9^ii\9)mini'^J2\9)m2n2^i3\9)m3n3 , , , 



Writing Trj[g ^)^^ = ■Kj{g*)nm = '^*ji9)nm in the second integral of (A.l), one gets the complex 
conjugate of the 3jm symbols defined above. Thus the contribution to the inner vertex is as 
mentioned in (4.11). 

B Tracing 

We briefly describe the tracing mechanism at the Inner vertex. To trace over the basis states 
of iJiiTnnnii) = \/2jji + ^T^ju{h)munu we simply use the inner product of {jmn\j'm'n') = 
Sjj'Smm'^nn'- Once the Kronecker delta functions have been implemented; the wavefunction at 
the inner vertex is of the form 

y^(27 +l)(-^^^ ^^ Jk\fjn JH 3k\ f jk 312 jn\f3k 312 i/sA 

3 3 

'^ n ^J^' i9)mnnn JJ ^JI^ {9)m'j^n'j^^n'j^nn^m'^^mi2^mi3m'j.j}PjH {9H)mHnHi'j'„ {9H)m'jjn'^ 
1=1 i=l 
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X 

3k' JI2 JI3\ fjk' jn JI3 



We then sum the 2>jm symbols with common entries, e.g. 



E( jk 312 3n\ ( 3k' 312 JI3\ _ hdk^n^ ■ ■ ■ \ 
U QI2 qi3 W 112 qi3 - {2h + l)^^'^''^''^' 

where {J1J2J3} is the 3j symbol and is such that it is 1 if j'l +J2 + J3 is an integer and \ji — J2I < 
J3 ^ ii + J2, and otherwise. This process then finds equation (4.12). 
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